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ABSTRACT. The first sources of ionizing radiation to condense out of the dark and neutral IGM sent 
ionization fronts sweeping outward through their surroundings, overtaking other primordial gas-clouds 
and photoevaporating them. Results are presented of the first gas dynamical simulations of this process, 
including radiative transfer, along with some observational diagnostics. 

1. Ionization Fronts in the IGM 



■ The role of hydrogen molecules as cooling agents in primordial gas clouds, necessary to 

form the first stars out of a gas of H and He with no heavier elements, is intimately 
related to the fate of those clouds in the presence of ionizing and dissociating radiation. 
That such radiation had an important effect on these primordial clouds is certain. The 
• neutral, opaque IGM out of which the first bound objects condensed was dramatically 

reheated and reionized at some time between a redshift z w 50 and z w 5 by the radiation 
Q ■ released by some of these objects (cf. [1] and references therein). When the first sources 

turned on, they ionized their surroundings by propagating weak, R-type ionization fronts 
which moved outward supersonically with respect to both the neutral gas ahead of and 
the ionized gas behind the front, racing ahead of the hydrodynamical response of the 
IGM [2,3]. The effect of density inhomogeneity on the rate of I-front propagation was 
previously described by a mean "clumping factor" q > 1, which slowed the I-fronts by 
increasing the average recombination rate per H atom inside clumps. This suffices to 
describe the rate of I-front propagation as long as the clumps are either not self-shielding 
or, if so, only absorb a fraction of the ionizing photons emitted by the central source. 
What is the dynamical effect of the I-front on the density inhomogeneity it encounters, 
however? 

The answer depends on the size and density of the clumps overtaken by the I-front. 
The fate of linear density fluctuations depends upon their Jeans number, Lj = A/Aj, 
the wavelength in units of the baryon Jeans length in the IGM at temperatures of 
order 10 4 K. Fluctuations with Lj < 1 find their growth halted and reversed (cf. [4]). 
For nonlinear density fluctuations, however, the answer is more complicated, depending 
upon at least three dimensionless parameters, their internal Jeans number, Lj = R c /Xj, 
the ratio of the cloud radius R c to the Jeans length A,/ inside the cloud at about 



10 4 K, their "Strom grcn number" L s = R c /i s , the ratio of the cloud radius R c to the 
Stromgren length £ s inside the cloud (the length of a column of gas within which the 
unshielded arrival rate of ionizing photons just balances the total recombination rate), 
and their optical depth to H ionizing photons at 13.6 eV, th, before ionization. If th < 1, 
the I-front sweeps across the cloud, leaving an ionized gas at higher pressure than its 
surroundings, and exits before any mass motion occurs in response, causing the cloud 
to blow apart. If th > 1 and L s > 1, however, the cloud shields itself against ionizing 
photons, trapping the I-front which enters the cloud, causing it to decelerate inside the 
cloud to the sound speed of the ionized gas before it can exit the other side, thereby 
transforming itself into a weak, D-type front preceded by a shock. Typically, the side 
facing the source expels a supersonic wind backwards towards the source, which shocks 
the IGM outside the cloud, while the remaining neutral cloud material is accelerated 
away from the source by the so-called "rocket effect" as the cloud photoevaporates 
(cf. [5]). As long as Lj < 1 (the case for gas bound to dark halos with virial velocity 
less than lOkms -1 ), this photoevaporation proceeds unimpeded by gravity. For halos 
with higher virial velocity, however, Lj > 1, and gravity competes with pressure forces. 
For a uniform gas of H density nu, c , located a distance tm P c (in Mpc) from a UV 
source emitting N p h^Q ionizing photons (in units of 10 56 s _1 ), the Stromgren length is 
only t a = (50 pc)(A r p h.56/?'Mpc)( n H.c/0.1 cm -3 ) -2 . We focus in what follows on the self- 
shielded case which traps the I-front. Some of these results were summarized previously 
by us in [6]. 

2. The Photoevaporation of an Intergalactic Cloud Overtaken by a Cosmo- 
logical Ionization Front 

As a first study of these important effects, we have simulated the photoevaporation of 
a uniform, spherical, neutral, intergalactic cloud of gas mass 1.5 x 1O 6 M , radius R c = 
0.5 kpc, density njj.c = 0.1 cm~ 3 and T = 100 K, in which self-gravity is unimportant, 
located IMpc from a quasar with emission spectrum F v oc v~ XA {v > ^h) and iV p h = 
10 56 s _1 , initially in pressure balance with an ambient IGM of density 0.001 cm -3 which 
at time t = has just been photoionized by the passage of the intergalactic R-type 
I-front generated when the quasar turned on. [A standard top-hat perturbation which 
collapses and virializes at z co ii = 9, for example, with total mass = 10 7 M Q , has circular 
velocity v c = Tkms" 1 , R c = 560 pc, and nn. c = 0.1 cm~ 3 , if Cl^h 2 = 0.03 and h = O.5.] 
Apart from H and He, the cloud also contains heavy elements at 10~ 3 times the solar 
abundance. Our simulations in 2D, axisymmetry use an Eulerian hydro code (called 
CORAL), with Adaptive Mesh Refinement and a Riemann solver based on the Van Leer 
flux-splitting algorithm, which solves nonequilibrium ionization rate equations (for H, 
He, C, N, O, Ne, and S) and includes an explicit treatment of radiative transfer by taking 
into account the bound- free opacity of H and He [7,8,9]. Our grid size in cylindrical 
coordinates (r, x) was 128 x 512 cells (fully refined). 

Most current simulations of cosmological gas dynamics which include photoioniza- 
tion do so by approximating the radiation field as uniform and isotropic, ignoring the 
inhomogeneity of the radiation field and setting optical depth to zero. To illustrate the 
importance of a more realistic approach to radiative transfer, we present results for two 




Fig. 1. ZERO OPTICAL DEPTH APPROXIMATION. One time-slice 50 Myr after turn-on of 
quasar located 1 Mpc away from cloud to the left of computational box along the a;-axis. From 
top to bottom: (a) isocontours of pressure, logarithmically spaced, in (r, x)— plane of cylindrical 
coordinates; (b) pressure along the r = symmetry axis; (c) temperature; (d) H I fraction; (e) 
He I (solid) and He II (dashed) fractions. 



cases in what follows: (1) zero optical depth and (2) optical depth properly included. 
Figure 1 shows the structure of the cloud 50 Myr after it was overtaken by the quasar's 
I-front as it sweeps past the cloud in the IGM when optical depth is neglected. The 
cloud and IGM are both instantaneously photoionized everywhere in this case, causing 
the overpressured cloud to expand isotropically into the surrounding IGM, acting as a 
spherical piston which sweeps up the IGM and drives a shock into it. The cloud matter 
eventually expands as a shell and evacuates a spherical hole. The column densities of 
H I, He I and II, and C IV for cloud gas of different velocities as seen along the symme- 
try axis at different times are shown in Figure 3. This cloud would initially resemble a 
velocity-broadened Lyman alpha forest quasar absorber ("LF") 10's of kms -1 wide, with 
A r Hi > 10 16 cm~ 2 , which evolves toward a narrower LF absorber with JVh,i < 10 15 cm~ 2 , 
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Fig. 2. Same as Figure 1, except: OPTICAL DEPTH INCLUDED. Bottom panel is bound-free 
optical depth along r = axis at the threshold ionization energies for H I (solid), He I (dashed), 
He II (dotted). 



with Nu e ii/Nm < 10 and New ~ 10 cm -2 throughout. As the spatial variations 
for the relative abundances of selected metal ions at 50 Myr plotted in Figure 5 show, 
only highly ionized metals are present in this case. By contrast, we show results for all 
these same quantities from the simulation which takes proper account of optical depth, 
instead, in Figures 2, 4, and 6. Since is <C R c initially, the cloud traps the I-front, as 
described above, and drives a supersonic wind from the surface facing the quasar. It 
takes more than 100 Myr to evaporate the cloud, accelerating it to 10's of kms -1 in the 
process. Key features of the flow are indicated by the numbers which label them on the 
temperature plot in Figure 2:1 = IGM shock; 2 = contact discontinuity between shocked 
cloud wind and swept up IGM; 3 = wind shock; between 3 and 4 = supersonic wind; 
4 = I-front; 5 = shock preceding I-front; 6 = shock that leads the motion of remain- 
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Fig. 3. Cloud column densities (cm 2 ) along symmetry axis at different velocities: ZERO OP- 
TICAL DEPTH APPROXIMATION. From top to bottom: (top) H I; (middle) He I (solid) 
and He II (dotted); (bottom) C IV. Each box labelled with time (in Myrs) since QSO turn-on. 
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Fig. 4. Same as Figure 3, except OPTICAL DEPTH INCLUDED. 




Fig. 5. Carbon, Nitrogen, and Oxygen Ionic Fractions Along Symmetry Axis at t = 50 Myr: 
ZERO OPTICAL DEPTH APPROXIMATION. 



ing neutral cloud gas into the shadow region. At early times, the cloud gas resembles a 
weak Damped Lyman Alpha ("DLA") absorber with small velocity width (<~ lOkms -1 ) 
and .ZVhi ~ 10 20 cm -2 , with velocity-broadened LF-like wings (width ~ 20kms~ 1 ) with 
A^hi ~ 10 16 cm~ 2 on the side moving toward the quasar, with a C IV feature with 
New ~ 10 12 cm~ 2 displaced in this same asymmetric way from the velocity of peak H I 
column density. After 160 Myr, however, only a narrow H I feature with LF-like column 
density N H i ~ 10 14 cm~ 2 remains, with N HcU /N H i ~ 10 2 and N C iv/Nm ~ [C]/[C] . 
Unlike the zero-optical-depth approximation, Figure 6 shows the presence at 50 Myrs 
of low as well as high ionization stages of the metals. 
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Fig. 6. Same as Figure 5, except OPTICAL DEPTH INCLUDED. 
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